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ABSTRACT 

We study the formation of fifty-three galaxy cluster-size dark matter halos (M = lO'"* ""'"^'^^Mq) formed within 
a pair of cosmological ACDM N-body simulations, and track the accretion histories of cluster subhalos with 
masses large enough to host ^ O.IL* galaxies. By associating subhalos with cluster galaxies, we find the 
majority of galaxies in clusters experience no "pre-processing" in the group environment prior to their accretion 
into the cluster. On average, ^^ 70% of cluster galaxies fall into the cluster potential directly from the field, 
with no luminous companions in their host halos at the time of accretion; and less than ~ 12% are accreted as 
members of groups with five or more galaxies. Moreover, we find that cluster galaxies are significantly less 
likely to have experienced a merger in the recent past (< 6 Gyr) than a field halo of the same mass. These 
results suggest that local, cluster processes like ram-pressure stripping, galaxy harassment, or strangulation 
play the dominant role in explaining the difference between cluster and field populations at a fixed stellar mass; 
and that pre-evolution or past merging in the group environment is of secondary importance for setting cluster 
galaxy properties for most clusters. The accretion times for z = cluster members are quite extended, with 
^ 20% incorporated into the cluster halo more than 7 Gyr ago and ^ 20% within the last 2 Gyr By comparing 
the observed morphological fractions in cluster and field populations, we estimate an approximate timescale 
for late-type to early-type transformation within the cluster environment to be ~ 6 Gyr 

Subject headings: cosmology: theory, large-scale structure of universe — galaxies: formation, evolution, high- 
redshift, interactions, statistics 
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One suggestion is that "pre-processing" in the group envi- 
ronment prior to cluster formation is import ant in setting the 
early-type fraction in clusters (e.g., Zablu doff & Mulchaevl 
119981; IZabludoffll2002l) . A related possibility is that galax- 
ies in clusters typically experience more mergers than field 
galaxies (prior to their accretion), and that this merger 
history bias plays a role in ex plaining population differ- 
ences dToomre & Toomr3ll972l) . Finally, the fact that the 
overall mix of galaxies in cluste rs by type is kno wn to 
evolve with reds hift (Butc her & Oem ler 1978; Ellingso n et alj 
SpoT; 'Tran et al. 2005; Gerke et al. 2007; Capak et al. '2007'; 
Coil et al. 2008; Loh et al. 2008; Finn et al. 2008) suggests 
that the internal cluster processes play a major role in set- 
ting the differences between the cluster and field populations. 



Here we examine the formation of clusters in a ACDM cos- 
mological simulation in order to gain insight into these ques- 
tions. 

Galaxies in clusters and groups are subject to a number of 
processes that may suppress star formation or change the mor- 
phology of a galaxy . One such effect, r am-pressure stripping 
(iGunn&Gotfl 19721; lOuilis et al J 120001) seems to have been 
observed directly in Virgo dChung et al J 120071) . and simula- 
tions suggest that this process should operate on short, ^ 1 
Gyr, timescales ( Tonnesen 2007) . Other processes of rele- 
vance includ e galaxy-galaxy "h arassment" within the clus- 
ter potential (iMoo re et al.lll99 6|) and cold g as "strangulation" 
dLarson et aljl9 80; Kauffmann et aljl993h . which cuts of the 
gas supply for ongoing star formation in cluster galaxies. 

By incorporating subsets of these expected cluster pro- 
cesses into semi-analytic plus N-body models, several 
groups have investigated cluster galaxy population trends 
in the context of AC DM (Baloghetal. 2000; Benson etaf] 
2001"; 'Diaferio et al. aOOlt lOkamoto & Na gashima ,200U 
Springel et al. 2001). Broadly speaking, these models have 
been successful in producing the general behavior that early 
type galaxies are more common in clusters, but full agree- 
ment between theory and observation has yet to be achieved. 
This is likely because some of the relevant processes (e.g., 
ram pressure stripping, harassment) have been neglected, and 
perhaps because some of the effects that were included (e.g., 
morphological transformation via mergers) were modeled by 
rough approximations. Our approach is related to these past 
theoretical efforts, but different in its goal. Specifically, we 
aim to quantify the basic assembly statistics for cluster mem- 
ber halos using N-body simulations and to present these as a 
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basis for evaluating scenarios and interpreting observations. 
Surely, many of the of the statistical results presented here 
were implicitly included in past theoretical models that used 
N-body simulations as a basis, but our aim is to present the 
ACDM predictions as purely as possible, without obscuring 
them with any particular set of model assumptions for the 
baryon physics. 

In what follows we use a pair of N-body simulations to 
track the assembly history of cluster-size dark matter halos. 
In section §|2]we discuss the simulations and the method for 
finding halos and subhalos. We present our findings in § [3] 
We reserve §|4]for a discussion of potential impUcations. We 
conclude in §|5] 

2. METHODS 

We study the formation histories of fifty-three M > 
1O''*/!~'M0 cluster-size dark matter halos extracted from two 
cosmological N-body simulations with comoving cubic vol- 
umes of 120/i~'Mpc and 80/;~'Mpc on a side. Each simulation 
corresponds to a flat ACDM cosmology with fl^ = 1 - ^a = 
0.3, h = 0.7, and cg = 0.9 and were performed using the Adap- 
tive Refinement Tree (ART) co d e of |Kravtso v et al. ( 199^. 
As discussed in lAUgood et alJ (l2006l) and .Wechsler et aP 
(12006), the 80/i~^Mpc simulation followed the evolution of 
5 12-' particles with a mass of 3.18 x lO^h'^M© and achieved 
a maximum force resolution of 1.2/i~'kpc. The 120/i~'Mpc 
simulation (lAllgood et al.l l2006h followed the evolution of 
5 12-' particles with a mass of 1.07 x 10^h"'MQ with a maxi- 
mum force resolution of 1.8/i~'kpc. 

We identify halos in the simulation using a varia- 
tion of the Bound Density Maxima Algorithm (BDM 
iKlvpin et al.l 1 1999b. spe cificallv adopting the methods out- 
lined in iKravtsov et al.l (12004.) . As described in Stewart et 
al. (2008) , virial radii and virial masses for halos are set by 
the radius within which the average density is Ay j, times the 
mean density of the universe (e.g.. iBrvan & Nor man 1998). 
At z = this definition implies a radius-mass relation of 
Ryi, ~ 951(Mvi,-/10''*/i-'Mo)'/l Subhalos are defined to be 
self-bound halos with centers located within the radius of a 
larger halo. In most cases, subhalo density profiles become 
overwhelmed by the larger halo's density field at a truncation 
radius, R,, that is smaller than /?vii- The truncation radius is 
defined to be the radius where a halo's density profile flattens 
to a value larger than -0.5. We define a halo's mass, M, to be 
the minimum Mvi, and the mass within Rf. 

We use 48 snapshot outputs from the 80/i~'Mpc simulation, 
spaced roughly equally in expansion factor a = 1/(1 +z) back 
to z = 2 1 .6 to generate merger histories for halos and subhalos. 
For the 120/i~'Mpc simulation we used 91 snapshot outputs, 
spaced roughly equally in expansion factor back to a z = 10. 1 . 
The merger trees were derived using the methods discussed 
in lAllgood et al.1 (l2006h . IWechsler et all (l2006h . and Stewart 
et al. (2008) . When halos are accreted into larger halos, we 
record the mass at this time and label it M,„. For halos that 
are not subhalos (field halos) we set M,„ = M. We use this 
mass, Mjn, as a proxy for luminosity in defining our cluster 
galaxy samples (see below). Once a halo becomes a subhalo, 
we continue to track its mass as it evolves within the larger 
halo. When subhalos fall below a critical mass, Mcr < Mj„, 
we explicitly remove them from our catalogs, and assume that 
any galaxy it was hosting has fallen out of the observational 
sample (either because it has lost a significant fraction of lu- 
minous mass or because it has been disrupted). This choice 
for Mcr allows us to define a sample cleanly at a mass scale 



where our halo finder is complete and our simulations are not 
strongly affected by over-merging. 

In what follows we will assume that halos and subhalos 
with Mj„ larger than a specific threshold will host one galaxy 
at their center. If a halo contains one or more subhalos, it is 
assigned one "central" galaxy in addition to one galaxy for 
each of its subhalos. The term host is used to describe the 
largest halo that contains a galaxy. A halo need not contain 
a subhalo in order for it to be classified as a host, it simply 
needs to be massive enough to contain a galaxy at its center. 
A subhalo cannot be a host. Finally, we term field halos to be 
all halos that are not contained within a larger halo. Thus, by 
definition field halos cannot be subhalos. 

By associating galaxies with subhalos larger than a critical 
mass at the time of their accretion, we are adopting a strat- 
egy similar to that used successfully by Conr oy et al.l (|2006|) 
and iBerrier et al.l (l2006l) . These authors were able to repro- 
duce both the large-scale and small-scale clustering statistics 
of galaxies by assuming a monotonic relationship between the 
luminosity of a galaxy and the maximum circular velocity that 
its halo had when it is first accreted into a larger halo (see 
IWang et al] |2006, for a similar approach). 

Our primary population of cluster galaxies is defined by set- 
ting Mjn > 10''-^ and Mcr > 1O''"/i"'M0. Averaging over 
both simulation volumes, this choice defines a sample with 
number density rig = 0.012/!-'Mpc~ . Matching this num- 
ber density to the the Slo an Digital Sky Surv ey (SDSS) r- 
band luminosity function (iBlanton et al.ll20()3h . we estimate 
that these cuts correspond to a galaxy population with an r- 
band magnitude brighter than M,- ~ -18.5 for h = 0.7. This 
is comparable to the lum inosity ranges used in most cluster 
morphology studies (e.g.. lDressledll980tlDressler et al.lll997t 
iPostman et al.ll2005h . Note that our results are insensitive to 
the precise choice of Mcr- We have redone the analysis de- 
scribed below and using an Mcr value that differs by a factor 
of ~ 3 from our fiducial Mcr and find virtually identical re- 
sults. 

Twelve of our cluster halos are taken from the higher reso- 
lution 80/i~'Mpc simulation. For these high-resolution clus- 
ters, we explore the accretion histories of a second set of 
lower-mass cluster subhalos with M,„ > lO"*^ and Mcr > 
1O"''^/i~'M0. This sample has a number density of rig = 
0.042 /i-^Mpc, which, with h = 0.7, corresponds to minimum 
lumin osity of Mr c^ -16.4 using the SDSS luminosity func- 
tion (IBlanton et alJl2003h . This sample will be used to test 
the dependence of our results on sample selection and should 
represent galaxies that are comparable to (although somewhat 
fainter than) the observational sample studied by iTreu et al] 
(2003). 

Of course, these estimated cluster galaxy luminosities are 
only approximate. We certainly don't expect that a simple 
mapping between mass and luminosity will hold in detail, 
especially within the cluster environment. However, as we 
show below, our results do not depend strongly on the adopted 
mass cut. This suggests that the uncertain mapping between 
individual halo masses their associated galaxy luminosities 
should not hinder the interpretation of our results. 

The total sample of 53 cluster halos have z = masses span- 
ning Mcius = (1.0-5.8) X IO^'^Ii'^Mq, with a median mass of 
1.48 X 1O''*/;~'M0. The total number of galaxy subhalos in 
this, our main sample of M,„ > lO"^/i~'M0 subhalos is 834 
(^ 15 per cluster). We also explore trends with cluster halo 
mass by dividing our main cluster sample into three mass bins 
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containing 18, 17, and 18 clusters each, with mass ranges that 
span 10i4"-i4.i3^ 101413-14.24^ ^nd 10i4-24-i4.76;j-iyi^g^ respec- 
tively. The clusters in these respective samples contain an av- 
erage of ~ 9, '--^ 14, and ^ 25 galaxies each. The subset of 12 
high-resolution cluster halos we study have z = masses that 
span Mdus = (l -0-3.4) X lO^'^h'^MQ. These clusters host a to- 
tal of 643 subhalos that meet ourM,„ > IQU'Vi'IMq criterion 
and have an average of ^ 54 galaxies per cluster. The most 
massive cluster in this sample hosts 102 galaxies. A summary 
of the clusters we study along with some properties of their 
accretion histories are given in Table 1 . 

In order to investigate the assembly of more massive 
clusters we also study cluste rs generated s emi-an alytically 
using the code described in IZentner et al.l (I2005L see also 
IZentner & BullocM l2003h . This model uses the extended 
Press-Schechter (EPS) formalism, applying the specific im- 
plementation of ISomerville & Kolatj (1 19991) to produce mass 
accretion histories for host halos of a given present-day mass, 
and then tracks the orbital evolution of each tree's subha- 
los in their host potential via analytic prescriptions govern- 
ing dynamical friction and tidal mass loss. The initial or- 
bital energy and impact parameter for each satellite system 
are drawn from probability distributions yielded by cosmo- 
logical N-body simulations. This algorithm produces sev- 
eral statistics which match those produced in well understood 
cosmological simulations. While not as accurate as our nu- 
merical simulations, this computationally-inexpensive semi- 
analytic model allows us to explore trends with cluster mass 
(from IO'^/j-IMq < M^us < lO^^-^h-^MQ) with thousands of 
stochastically-generated halo formation histories. 

3. RESULTS 
3.1. The environmental history of cluster galaxies 

Before presenting results on the characteristics of cluster 
galaxies, we first explore the mass assembly distribution for 
the clusters themselves. The solid, monotonically increasing 
line in Figure[T]shows the mass fraction of our z = sample of 
53 clusters that was built up from accreting halos more mas- 
sive than M, where M is the halo mass just prior to accretion 
into the cluster Note that we have normalized the accreted 
halo mass M relative to the final cluster mass at z = 0, M/McUs- 
The histogram shows the corresponding differential distribu- 
tion. We see that the distribution peaks at M/M^us ~ 0.1. 
This result is consistent with the expectation that dark matter 
halos of any size mass are built primarily from objects that 
are ~ 10% o f the mass of the fina l halo (Purcell et al. 2007; 
[Stewart et aLll2008l: IZentneill2007l) . Note that more than half 
of the mass accreted from objects large enough to host galax- 
ies (M/Mcius ^ 0.001 Mcius) is accreted from group-scale sys- 
tems with M/Mcius > O.lMcius- 

Compare this result to Figure |2] which shows the fraction 
of z = cluster galaxies that were accreted as members of host 
halos of a given mass. The short-dashed (blue), solid (black) 
and long-dashed (red) lines correspond to clusters in our three 
mass bins, centered on z = masses Mdus — lO''*"^, lOi'*'^, and 
1014.35^-1^^^ respectively. We see that for the typical cluster 
in our sample, only ^25% of the cluster's galaxies were ac- 
creted as part of group-size objects with M > 10^^ Ii'^Mq. For 
the lowest and highest mass subsamples we see that ~ 15% 
and ~ 30% of galaxies are accreted from group-mass halos. 
The results are not sensitive to the selection of our M^.,- value. 

The previous two figures show that cluster mass assembly is 
dominated by the most massive (group-size) accretion events, 
while cluster galaxy assembly is dominated by lower-mass 
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Fig. 1. — Cumulative and differential fraction of cluster mass, Mdus, ac- 
creted in halos of mass M/M^ius over the history of the simulation. This is 
an average result based on our sample of 53 clusters with a typical mass of 

Mdus~ 10'4.2 /r'M0. 



SO.I 



fe 0.6 



,0.4 



t. 0.2 



M„,„. =^10"'-"is h-' M„ 




1012 io'3 lO'" 
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Fig. 2. — Cumulative and diiferential fraction cluster galaxies that fell 
into their respective clusters as part of a host halo of a given mass. Cluster 
galaxies were identified with M,„ > 1O"'^/i"'M0, con'esponding to galax- 
ies brighter then Mr ~ —18.5. The red-dash, black-solid, and blue-dot-dash 
lines correspond to increasing cluster halo mass bins, as labeled. The lines 
that increase monotonically towards lower accreted halo masses are the same 
distributions presented cumulatively. 



(galaxy-size) halos. This can be understood by noting that the 
number of galaxies that a halo hosts does not increase linearly 
with host mass. A small group (M ^ IQ^^'^h'^MQ) contains ~ 
10 times the mass of a single-galaxy halo (M ^ 10^^'^ Ii'^Mq). 
However, typically a group halo of this mass will host only ~ 
2-3 galaxies that are as b right as the ga laxies associated with 
1O11'^/!~1M0 halos. (e.g.. iBerrier et alJ l2006. and references 
therein). This means that a small number of group-size halos 
can deposit a significant amount of mass into a cluster without 
contributing an equally large fraction of its galaxies. 

While the mass of a galaxy's host at accretion provides 
some insight into the environment within which it formed and 
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Fig. 3. — The fraction of surviving cluster galaxies, A'g/A'totaii that fell 
into their clusters along with a given number of companions in their host 
dark matter halo. The black solid line corresponds to our standard galaxy 
sample (M,-^ - 18.5) while the green dashed line shows the same distribu- 
tion for our high resolution sample (M,„ > 1O"/?"'A/0, M, < - 16.4). The 
inset box shows the distribution out to the maximum number of companions 
found in our sample in order to illustrate the tail of the distribution. Note that 
"companions" do not have to survive in order to be included in the count. 
We see that, on average, r^ 70% of cluster galaxies fell into their clusters 
directly from the "field" (with zero companions in their dark matter halos), 
while ~ 80% fell in with 2 or fewer companions. 



evolved before joining the cluster, a more direct measure of 
the environment can be obtained by counting the number of 
galaxies that each host halo contained at the time of its infall. 
Figure [3] shows the distribution of the number of companions 
each cluster galaxy had within its host halo at the time it was 
accreted into the cluster The companion count includes any 
other galaxies that existed within a galaxy's host halo at the 
time of accretion. We do not require that a companion galaxy 
"survives" with M > M^r at z = in order to include it in this 
count — it simply must have M > Mcr at the time of accretion 
into the cluster A companion count of zero implies that the 
galaxy was accreted as the only object in its halo {i.e., it was 
accreted from the field). The two solid lines show the binned 
and cumulative distributions for our standard sample of 834 
galaxies within 53 clusters, and the pair of dark green dashed 
lines correspond to our sample of 643 lower mass galaxy ha- 
los taken from the 12 high-resolution clusters. For our stan- 
dard sample, approximately 70% of cluster galaxies were ac- 
creted as the only galaxy within their halo and ^ 88% were 
accreted with fewer than 5 galaxies. The numbers are similar 
for the sample that includes smaller galaxy halos. In this case, 
^ 65% of the cluster galaxies are accreted directly from the 
field. 

So far we have considered only averages for entire sam- 
ples and sub-samples of our cluster-galaxy population. We 
would also like to obtain some indication of the variation in 
assembly histories from cluster to cluster. Figure H] presents 
statistics on the individual (M,„ > 10"'^/i~'Mo) galaxy pop- 
ulations subdivided into their respective clusters. The red- 
dashed line shows the distribution of our clusters that have 
a given fraction (/) of their galaxies accreted directly from 
the "field", as the only object in their halos. For example, 
an abscissa value of / = 0.5 implies that 50 % of a cluster's 
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Fig. 4. — Distribution of galaxy clusters with agiven assembly history. The 
red-dashed line shows the frequency of simulated clusters that have a given 
fraction / of their galaxies accreted directly from the field (with no compan- 
ions in their host halos as they fell in, / = A'o/Mo(al)- The red-dashed peak at 
/ = 1 shows that ^^ 18% of our clusters have had 100% of their galaxies ac- 
creted directly from the field. The blue-dot-dashed line shows the frequency 
of clusters with a given fraction of their galaxies accreted with five or more 
companions in their host halos (/ = A'>5 /Wtotal)- The large peak at / = shows 
that ~ 83% of our simulated clusters have had none of their galaxies accreted 
from groups with five or more companions. Finally, the solid line shows the 
frequency of simulated clusters with a given fraction of their galaxies that 
were accreted as part of a group-mass halo with M > IO'^/i^'Mq. The broad 
bump around / ~ 0.2 shows that, typically, our clusters have ~ 10—40% of 
their galaxies accreted from halos larger than the ~ IO'^/i^'Mq. Here we 
used our standard galaxy sample with M,„ > lO"-^/r'M0, M, < - 18.5. 



galaxies were accreted from the field. The peak in the red- 
dashed histogram at / = 1 shows that ~ 18% (10/53) of our 
clusters were assembled entirely from field accretions. The 
blue, short-dashed line shows the distribution of our clusters 
that had a fraction / of their galaxies accreted together with 
5 or more companions in their host halos at the time of their 
infall into the cluster The spike in the blue short-dashed his- 
togram at / = implies that ~ 83% (44/53) of our clusters 
had none of their galaxies accreted from a group with 5 or 
more companions. Finally, the solid line shows the frequency 
of clusters with a fraction / of their galaxies accreted from 
within group-mass halos with M > lO'-^/z'^M©. We see that 
most clusters have between / = 10% and 40% of their galax- 
ies accreted from group-mass halos. However, the solid-line's 
peak at / = shows that a non-negligible fraction (^ 28%, 
15/53) of our clusters had none of their galaxies accreted as 
members of group or cluster-size halos. 

3.2. Trends with Cluster Mass 

Figure ID shows that more massive clusters are more likely 
to have had their galaxies accreted from massive host halos. 
Figure |5] explores the issue of mass dependence more fully. 
Shown is the fraction of surviving cluster galaxies that fell 
into the cluster without a bright companion in their host halo 
at the time of accretion as a function of the cluster mass. The 
three triangle data points show the median fractions for three 
mass bins based on our standard sample of 53 clusters with 
galaxy halos set via M,„ > 10"-^/!"'Mo, Mr< - 18.5 . The 
square data point is from our high-resolution sample, with 
galaxies identified using ourM,„ > 10" /!~'Mq criterion, cor- 
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responding to M,-;$ - 16.4. Vertical error bars correspond to 
the sixty-eighth percentile range and the horizontal error bars 
span the mass range of clusters included in the bins. As ex- 
pected, we do see a slight mass trend in the median, with 
the fraction of field-accreted galaxies falling from ~ 75% to 
^ 65% as the cluster mass increases from ^ 10'"^/!~'Mq to 
~ 3 X IQ^'^H'^Mq, although the overall trend is weak com- 
pared to the scatter from cluster-to-cluster at fixed mass. 

In order to explore how this trend con tinues over a broader 
cluster mass range, we have used the Zentner et al.l ( l2005h 
semi-analytic mer ger-tree and substructure code (see also, 
iPurcell et al.ll2007l) to generate 1000 cluster halo realizations 
at each of 10 cluster halo mass values between M^us = 10'^ 
and 2 x lO'''/i~'M0. The three line types correspond to 
three different choices for defining the galaxy samples, with 
M,„ > 10'^, 10"-^, and 10"/!"'Mq, from top-to bottom, re- 
spectively. Galaxy subhalos are counted at z = as surviving 
galaxies if they maintain masses larger than Mcr = 10", 10", 
and 1O"^'^/i~'M0, respectively. The long-dashed red and 
solid black lines should correspond to the red triangle and 
black square N-body points, respectively. While the overall 
normalization of the semi-analytic estimates is low compared 
to the N-body result, the mass trend is in approximate 
agreement, with more massive clusters accreting a smaller 
fraction of their galaxies directly from the field, with a mass 
scaling as ^ M'^l^^ over the mass range plotted. Note that 
we expect even the most massive clusters to have a signif- 
icant fraction of their galaxies accreted directly from the field. 



3.3. Cluster assembly with time 

The results presented in the previous section suggest that 
most cluster galaxies experienced no evolution in a group en- 
vironment prior to their accretion into the cluster, lending sup- 
port to the idea that internal cluster processes are responsible 
for the differences seen between cluster galaxies and those in 
field environments. If this is so, then the distribution of time 
spent in the cluster environment can provide insight into the 
timescales required for morphological transformation or the 
truncation of star formation. 

Figure |6] shows the distribution of accretion times for sur- 
viving cluster galaxies in our main cluster sample. We see that 
the accretion rate has been fairly uniform over the past ^ 8 
Gyr, with the median lookback time to accretion at ^^ 4 - 5 
Gyr Interestingly, as shown by the red-dashed line, galax- 
ies that were acccreted as part of group-mass halos with 
M > 1O'^/!~'M0 are biased to having been accreted later than 
the full sample. Specifically, the median lookback time for 
galaxies accreted as part of a group-mass system is ^^ 3 - 4 
Gyr This is not surprising, since it takes time for group-mass 
systems to form in a hierarchical universe. Indeed, we suspect 
that this fact is driving the double-peaked signature in the ac- 
cretion time histogram, which has a slight dip at ^ 4 Gyr. We 
do not find a trend between the accretion time distribution and 
cluster mass within our sample. Table 1 lists the median look- 
back accretion times for each cluster individually. 

Another way to gain insight into the history of cluster galax- 
ies is to quantify the amount of time a galaxy spends within 
a host halo of a given mass. Figure |7] shows the mean time 
that a cluster galaxy has spent in a host halo of a given mass, 
averaged over the history of the simulation. We see that this 
time-weighted host mass distribution for cluster galaxies is 
bimodal - on average, cluster galaxies have spent time ei- 
ther in their cluster or within a galaxy-mass halo prior to 
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Fig. 5. — Fraction of surviving cluster galaxies that fell into the clus- 
ter without a bright companion in their host halo at the time of accretion. 
The three triangle data points are taken from our main sample of simu- 
lated, where galaxy halos are defined withM;,, > 10"-^/i"'Mq, M,- < -18.5. 
The square data point is from our high-resolution sample, M,- < —16.4, 
with galaxies identified using our M,„ > 1O"/i"'M0 criterion. En'or bars 
along the y axis for the points reflect the 68 percentile range, and the er- 
ror bars along the x-axis reflect the mass range of clusters included in the 
bin. The line s correspond t o data generated by the semi-analytic substruc- 
ture code of Zentner et alj 420051) . The black (solid) line is for a sample 



of galaxies with M,„ > lO""/i''M0 (Mr. 



16.4), the red (dashed) line 
18.5), and the 



is for our standard mass cuts of M,,, > 10"^/i 'Mq (Mr^ — 
blue (dot-dashed) line is for a sample with M,>, > IO'^^/i^'Mq. Both the 
M„, > 10"5andM„, > lO'^" samples have a criteria of M„ > lO""/r'M0. 
The M,„ > lO" " sample matches our high resolution sample with Mcr > 
2 X IO'^^/i^'Mo. Note that all three samples show the same trends. The 
only difference between them is that the trend is simply offset based on the 
mass of the sample. 
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Fig. 6. — Differential and cumulative distributions of accretion times for 
our primary sample, M, <— 18.5 ofN-body cluster galaxies atz = 0. The solid 
black line shows the distribution of the whole sample. The red long-dashed 
line shows the cumulative accretion times for the subset of cluster galaxies 
that were accreted from group-mass halos with M > 10"/i"'A/q. Note that 
galaxies accreted from group halos ai'e biased to fall into clusters later than 
other galaxies. 
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accretion. This distribution is calculated by examining the 
mass of each galaxy's dark matter halo at each of the snap- 
shots taken from the simulation. The simplifying assump- 
tion is made that the halo spends all of its time at that given 
mass until the next timestep. The pronounced dip in the 
middle shows that cluster galaxies tend to avoid spending 
time in groups. Here we have used the full sample of clus- 
ter galaxies with M,„ > lO"'^/i~'M0 to compute the average 
time that a cluster galaxy has spent in a host halo of a given 
mass. The peaks of the distribution are at the single galaxy- 
halo scale, ~ 1O'''^/z~'M0, and at the median cluster scale 
^ 1O''*'^/i~'M0. Note that galaxies can spend some time in 
halos smaller than our M,„ threshold. Specifically, these halos 
grew from a small mass to a mass larger than 10"'^ Ii'^Mq be- 
fore being accreted. It is clear that, on average, cluster galax- 
ies spend very little time in the group-mass halos between 
lO'"andlO'^-^/i-'M0. 

Figure[8]makes a similar point, but in a more extreme fash- 
ion. Here we plot the fraction of galaxies that have spent any 
time within a halo within a given mass bin. We use the same 
simplifying assumption that the halo remains at the same mass 
until the next time step, therefore our ability to determine 
whether halos have spent anytime in a halo is limited by our 
output timestep spacing. In this figure, a value of unity im- 
plies that every galaxy has spent at least some time in a halo 
of this size. 

3.4. Cluster Galaxies and Merger Histories 

It is common to argue that spheroidal galaxies are associ- 
ated with halos that have undergone significant mergers. This 
motivates us to compare the merger histories of field and clus- 
ter galaxies. The left panel of Figure |9] shows the cumula- 
tive fraction of cluster galaxies halos (M,„ > lO"'^/i~'M0, 
Mr"^ - 18.5) that have have had a large merger since a given 
lookback time. The right panel shows the same statistic, now 
for^eW halos with M > lO"^/i"'M0. The different line types 
correspond to different mass-ratio mergers: m/Mo = 1/10 to 
2/5 from top (solid) to bottom (dashed). The mass ratio is 
defined to be the ratio relative to Mq = Mi„ for cluster galaxies 



and relative to Mq = M at z = for field halos. 

At large lookback times cluster galaxies and field galaxies 
show similar results, at least for intermediate-size mergers. 
The total fraction of systems with > 1/10 and > 1/5 mergers 
in the last ~ 10 Gyr are ^75% and ~ 40%, respectively, for 
both the cluster and field populations. The cluster population 
has a larger fraction of systems that have ever experienced 
very large mergers. For example, ~ 25% of cluster galaxies 
have experienced something larger than a 0.3 merger in the 
last 10 Gyr, compared to just ^ 15% in the field. 

Perhaps the most striking difference is that the cluster pop- 
ulation is much less likely to have had a recent (< 6 Gyr) 
merger event than galaxies in the field. The fraction of cluster 
galaxies with a significant merger in the last ^ 6 Gyr is less 
than ^ 5%. Given that half of cluster galaxies were accreted 
more than ^4-5 Gyr ago (Figure|6l), this general result is not 
surprising. We expect the high-speed cluster environment to 
greatly reduce the likelihood for a merger. Large mergers are 
more likely to occur in the field. 

4. DISCUSSION 
4.1. Pre-processing 

Galaxy groups provide an important intermediate en- 
vironment between the field and high-density clusters 
for testing ideas about galaxy formation and evolution 
(jZabludoff & Mulchaey 1998). Indeed, it is possible that pre- 
processing in the group environment prior to the assembly of 
galaxy clusters is an important factor in explaining why clus- 
ter galaxies differ significantly from the field population. As- 
suming that ACDM provides an accurate description of the 
universe, the simulations presented here allow us to character- 
ize the importance of groups in the global formation of galaxy 
clusters. 

Note that with our definition, groups and clusters are de- 
fined using the standard "virial" over-density boundary for 
host dark matter halos. In this case, most galaxies in the uni- 
verse do nof reside in groups. ACDM simulations suggest that 
only ~ 10% of ^L* galaxy hal os reside within the virial radii 
of group or cluster halos (e.g.. iBerrier et alj|2006l and refer- 
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ences therein). It is interesting then to ask whether majority of 
cluster galaxies were assembled from the field population or 
whether they are biased to be galaxies that evolved in groups. 

While Figure [T] demonstrates that a significant fraction of 
the mass in clusters is accreted from group-size dark matter 
halos, Figures |2] and [3] show that the same is not true for 
galaxies in clusters. Less than 17% of cluster galaxies in 
^ 1O'"^^/i~'M0 clusters fell in as a part of a group with five 
or more other galaxies. Similarly, ^ 25% were accreted as 
part of a halo more massive than 1O'^/;~'M0. This finding 
suggests that pre-processing cannot play the dominant role in 
differentiating cluster galaxy populations from the field. As 
we discuss below, our results suggest that global differences 
between cluster and field populations must be set by environ- 
mental influences associated with the clusters themselves. 

It is also important to emphasize that the pre-cluster evo- 
lution histories of cluster galaxies varies significantly from 
cluster to cluster. This variation is demonstrated in Figure |4] 
For example, while most of our clusters (44/53 ~ 83%) have 
galaxy populations that are completely devoid of objects that 
were accreted from groups with five or more companions, a 
small fraction (4/53 ~ 7.5%) of our clusters have a majority 
(/ > 0.5) of their galaxies that fell into the cluster from groups 
of this kind; 

Zabludof Tife Mulchaeyl (Il998l) have suggested that cD 
galaxies in clusters may form originally in the group envi- 
ronment. Our results do not prevent such a scenario from 
occurring in some cases, particularly in those where a large 
group-scale merger has occurred recently. As seen in Figures 
|7]and[8] there are always a few galaxies that evolve in groups. 
Additionally, these figures only include galaxies that survive 
to the present. Galaxies that are accreted into the cluster and 
destroyed may provide material for the growth of cD ga laxies 
and is consistant with the results of iLin & Mohrl (|2004|) . 

We also find a weak trend between the fraction of clus- 
ter galaxies that could have experienced pre-processing and 
the mass of the cluster. For the lowest mass third of our 



cluster sample, M^us — 10''*''^/!~'Mq, we see from Figure 
|5]that ^75% of cluster galaxies are accreted directly from 
the field. This number drops to ^65% in our largest mass 
sample, Mdus — lO''*-'^/i~'M0. This trend with mass is cer- 
tainly real, however, the variation from cluster-to-cluster is 
much stronger than t he mass trend itself. The results of 
IWeinmann et al.l (1200 6) suggest that the early-type fraction 
of cluster galaxies rises from ^ 50% to ~ 55% in clusters of 
mass ~ 10'^ to lO'^'^/z~'M0. It is possible that pre-processing 
could play a role in driving this weak mass trend. 

4.2. Cluster Galaxy Merger Histories 

As we showed in Figure |9] there are clear statistical dif- 
ferences between the merger histories of cluster galaxies and 
field galaxy halos. Nonetheless, it is unlikely that merger his- 
tories alone can explain the morphology, color, and spectral- 
type differences seen between cluster galaxies and the field. 
For example, let us take 60% and 30% as characteristic 
bulge-dominated fractions in clusters (within < /?vii) and the 
field, respectively (Whitmore & Gilmore 1991; Postman et al] 
l2005h . By examining Figure |9] we see that the field spheroid 
population potentially could be associated directly with the 
~ 30% of field halos that have experienced a > 1 /5 merger 
in the last ^ 9 Gyr. However, the 60% cluster spheroid popu- 
lation cannot be explained with the same set of assumptions: 
only '^ 25% of cluster halos have had a similar merger in the 
same period of time, and only ~ 45% have ever experienced 
a > 1 /5 merger. 

Given that merger histories alone cannot explain the dif- 
ferences between cluster and field populations, we are forced 
to construct more complicated scenarios involving mergers to 
help explain the observed population differences. For exam- 
ple, galaxies in the field are expected to be surrounded by 
reservoirs of baryons that can cool to reform a disk after a 
large merger Galaxies in clusters would likely be stripped of 
this reservoir by the ambient cluster medium, making it im- 
possible for cluster galaxies to accrete material to reform a 
disk after they have fallen into the cluster The right panel of 
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Figure|9]shows that we could potentially associate the ^ 30% 
of bulge-dominated galaxies in the field with the ~ 1/10 
mergers that occurred within the last ~ 6 Gyr (and assume 
that merger remnants, older than 6 Gyr, in the field have re- 
formed disks). However, even this scenario would have dif- 
ficulties explaining the cluster population: only ^ 70% of 
cluster galaxies have ever experienced a > 1/10 merger (left 
panel. Figure |9]l and only ^35% were accreted into the clus- 
ter in the last 6 Gyr Thus, obtaining more than a ^ 25% 
spheroid fraction without other transformational effects (i.e., 
harassment) would seem difficult. 

The scenario discussed above is si milar to that adopted 
in many s emi-analytic models (e.g., iBenson et al.l 120011: 
ISpringelet al. 2001), where galaxies experience "strangula- 
tion" as their fresh gas supply is cut off when they fall into 
the cluster environment. Not only does this effect prevent 
the possible reformation of disks, but it can cause the galaxy 
to redden over ^1-3 Gyr timescales tiPoggianti et al.ll 19991 : 
lBaloghetalJl2000t lEningsonetalJl2001h and thus help ex- 
plain the morphological and spectral mix of cluster popula- 
tions. It should be noted, however, that observations of mod- 
erate redshift clusters suggest that either two different physi- 
cal processes or at least two different timescales are responsi- 
ble for spectral and m orphological transformations in clusters 
(lPoggiantietalJll999l) . 

4.3. Accretion Times and Galaxy Transformation 

The timescales available for cluster galaxy transformation 
can be garnered directly from Figure |6] (also Figure |7]i. We 
expect gravitational processes like harassment to operate over 
timescales similar to the cluster dynamical time, l/^/Gpvir — 
4 Gyr. Given this, Figure|6]shows that ~ 70% of cluster galax- 
ies have been in the cluster potential long enough to experi- 
ence significant dynamical perturbations. In order to trans- 
form a ^ 70% field disk population to a ^^ 40% cluster disk 
population, we need ^^ 3/7 ~ 40% of accreted galaxies to 
be affected (assuming spheroids remain spheroids). Figure 
|6] shows that ~ 40% of cluster galaxies have been within the 
cluster longer than ~ 6 Gyr {^ 1.5 dynamical times). This 
estimate suggests that ~ 6 Gyr is a typical transformation 
timescale for cluster galaxies. 

The timescale for strangulation is much shorter than a dy- 
namical time, and should begin to operate as soon as the 
galaxy encounters the intra-cluster medium (< 1 Gyr). We 
expect ram pressure to be important at higher hot gas den- 
sities than strangulation (i.e., at smaller cluster-centric radii) 
and therefore the associated timescale for ram pressure to act 
should be somewhat longer than that for strangulation, per- 
haps of order a dynamical time. Given the infall time distri- 
bution in Figure |6] we expect that most cluster galaxies will 
have been affected to some degree by the cutoff of fresh gas 
infall (e.g., ^ 90% with infall times > 1 Gyr). This process 
can alter spectral properties without affecting the morpholog- 
ical mix. A smaller fraction have been in the cluster for ~ 1 
dynamical time {^ 50% with infall times > 4 Gyr) - long 
enough to experience significant dynamical (morphological) 
transformations. 

The above discussions provide some qualitative evaluation 
of ideas that may explain why the cluster galaxy population 
is different from the field population. A similar analysis, 
involving a more precise characterization of radial depen- 
dencies and merger histories, and the evolution of cluster 
galaxy infall times with the redshift of the cluster, will be an 
important avenue for future investigation. 



5. CONCLUSIONS 

We use two cosmological N-body simulations set within 
the concordance ACDM cosmology to study the formation of 
cluster-sized dark matter halos with masses spanning M^us = 
]^q14.o-i4.76^-i^^ Our primary results are based on track- 
ing the merger histories of galaxy subhalos within these clus- 
ters. These, our cluster galaxies, are picked to have masses 
M,„ > 1O"^/!"'M0 at the time they first fell into their host ha- 
los, and should correspond approximately to '^ 0.1 L* galax- 
ies. Our conclusions may be summarized as follows. 

• We find that the majority of cluster galaxies {^ 70%) 
are accreted directly from the field, as the only objects 
within their dark matter halos at the time of their infall 
into the cluster virial radius. A minority {^ 25%) were 
accreted as part of M > lO'-' /i~'M0 group-mass halos, 
and a small fraction (~ 12%) of cluster galaxies fell into 
their clusters with more than 5 galaxies in their halos at 
the time of accretion. 

• Cluster galaxy infall histories show significant variation 
from cluster to cluster. For example, 9/53 (~ 17%) of 
the clusters in our sample were assembled entirely from 
the accretion of field galaxies. However, 4/53 (~ 7%) 
of our clusters had a majority of their galaxies (/ > 0.5) 
accreted from groups with five or more companions. 
Therefore, while on average cluster galaxies tend to be 
accreted from the field, there are some clusters that do 
not follow this trend. 

• More massive clusters have a smaller fraction of their 
galaxies (at fixed luminosity) accreted directly from 
the field. Approximately 75% of ^ O.IL* galaxies in 
M ~ 1O'"^/!~'M0 halos are accreted from the field, com- 
pared to ^ 65% of galaxies in M ~ IQ^'^-^^H'^Mq clus- 
ters. The scatter in formation histories from cluster-to- 
cluster at fixed mass is more significant than this mass 
trend. 

• The median lookback time to accretion for galaxies 
within clusters is ^ 4.5 Gyr, and ~ 85% of galaxies are 
accreted between 1 and 9 Gyr ago. By assuming that 
cluster galaxies are accreted with a morphological mix 
similar to the field, we estimate an approximate cluster- 
environmental transformation timescale of ^^ 6 Gyr. 

• Galaxy subhalos in clusters are significantly less likely 
to have had a recent (< 6 Gyr) merger than similar mass 
galaxy halos in the field. The merger fraction within the 
past ~ 12 Gyr is ~ 5% higher for the cluster subhalo 
population. 

Taken together, these results suggest that the observed pop- 
ulation differences between galaxies in clusters and those in 
the field are driven primarily by internal cluster processes. 
Given ACDM as a basis, merging in the group environment, 
or any other type of pre-processing in galaxy groups prior to 
cluster assembly, cannot be a major factor in setting the nearly 
two-to-one difference in early-type fraction between clusters 
and the field. 

Approximately half of an average cluster's population is ac- 
creted more than 4 Gyr ago (~ 1 dynamical time). This allows 
ample time for gravitational processes to drive morphological 
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transformations within the cluster environment. Moreover, in- 
teractions with the intra-cluster medium that remove gas and 
suppress star formation in cluster galaxies likely begin to op- 
erate on even shorter timescales (^ 1 Gyr). Therefore, most 
cluster galaxies (~ 90%) will be affected by the cluster envi- 
ronment, at least to some degree. 

While our results suggest that pre-processing is not the 
dominant mechanism in setting galaxy cluster processes, we 
do find most clusters have a non-negligible fraction of their 
galaxies accreted from a group environment. This is espe- 
cially true for the more massive clusters in our sample (~ 
IO^'^-^^H'^Mq), for which we find that -- 28% of their galaxies 
were accreted as part of a halo larger than a group-mass scale 
M > 1O'^/!~'M0. Therefore, some amount of preprocessing 
should occur. 

As mentioned in the introduction, important constraints 
on the types of processes that act to shape the cluster and 
group galaxy populations come from studies at intermediate 
to high red shift (Butcher & Oemler 1978; Tranetal. 2005; 
iGerke et alJ l2007t iCapak et all 120071: iFinn et all 120081) . A 



valuable direction of future work will be to combine the pre- 
dicted N-body statistics for cluster halo assembly time with 
these results as a means to constrain specific scenarios for 
galaxy transformation and star formation suppression. 
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Table 1 
Cluster Halo Properties 



Cluster id 


Mdus 


Svir 


A'grp/Motal 


No/N,Ml 


N>5/N,aal 


Tin Mass 


Tj M Galaxies 


A'grp/A'total 




(lo'^r'A/e) 


(r'Mpc) 


(>10"5) 


(> 10"5) 


(> 10"'5) 


Gyr 


Gyr 


(> lO'iO) 


L12() 


5.86 


1.71 


15/53 


34/53 


6/53 


7.1 


2.4 


NA 


2.120 


4.41 


1.56 


10/31 


23/31 


0/31 


7.4 


5.0 


NA 


3.120 


3.80 


1.48 


3/23 


19/23 


0/23 


7.4 


4.9 


NA 


1.80 


3.44 


1.44 


6/28 


17/28 


0/28 


6.2 


5.2 


26/102 


4.120 


2.95 


1.36 


9/23 


18/23 


0/23 


4.9 


2.2 


NA 


5.120 


2.77 


1.34 


5/26 


14/26 


0/26 


5.2 


5.8 


NA 


6.120 


2.60 


1.31 


7/33 


21/33 


0/33 


5.2 


5.5 


NA 


7.120 


2.57 


1.30 


7/21 


14/21 


0/21 


3.4 


1.8 


NA 


2.80 


2.44 


1.28 


0/19 


19/19 


0/19 


7.8 


6.5 


6/89 


3.80 


2.11 


1.22 


3/19 


17/19 


0/19 


6.2 


3.0 


11/51 


8.120 


2.06 


1.21 


4/23 


9/23 


10/23 


8.4 


6.1 


NA 


9.120 


1.94 


1.19 


23/41 


21/41 


9/41 


6.6 


4.5 


NA 


10.120 


1.92 


1.18 


8/27 


16/27 


0/27 


4.7 


3.0 


NA 


11.120 


1.88 


1.17 


6/22 


18/22 


0/22 


8.1 


5.8 


NA 


4.80 


1.87 


1.17 


6/19 


13/19 


5/19 


5.8 


4.9 


25/71 


12.120 


1.85 


1.17 


2/19 


17/19 


0/19 


6.0 


5.5 


NA 


13.120 


1.80 


1.16 


2/17 


15/17 


0/17 


5.8 


4.5 


NA 


14.120 


1.74 


1.14 


12/16 


9/16 


7/16 


0.7 


1.3 


NA 


15.120 


1.71 


1.14 


0/17 


17/17 


0/17 


8.7 


5.7 


NA 


16.120 


1.60 


1.11 


6/16 


9/16 


0/16 


7.4 


2.1 


NA 


17.120 


1.59 


1.11 


7/15 


8/15 


0/15 


2.4 


3.0 


NA 


5.80 


1.55 


1.10 


0/18 


14/18 


0/18 


7.4 


3.0 


3/43 


6.80 


1.51 


1.09 


3/15 


10/15 


5/15 


6.2 


3.0 


9/47 


18.120 


1.51 


1.09 


5/13 


9/13 


0/13 


2.2 


2.2 


NA 


19.120 


1.49 


1.09 


0/10 


10/10 


0/10 


6.5 


7.1 


NA 


7.80 


1.48 


1.08 


6/23 


16/23 


0/23 


6.2 


3.9 


14/67 


20.120 


1.48 


1.08 


1/11 


6/11 


0/11 


5.8 


4.2 


NA 


21.120 


1.48 


1.08 


3/8 


7/8 


0/8 


3.6 


1.8 


NA 


22.120 


1.48 


1.08 


1/6 


4/6 


0/6 


7.6 


1.9 


NA 


23.120 


1.45 


1.08 


1/10 


8/10 


0/10 


7.3 


2.1 


NA 


8.80 


1.44 


1.07 


3/10 


5/10 


0/10 


6.8 


5.2 


9/42 


9.80 


1.42 


1.07 


1/9 


8/9 


0/9 


5.5 


3.9 


9/37 


24.120 


1.41 


1.07 


0/11 


9/11 


0/11 


7.6 


4.9 


NA 


25.120 


1.39 


1.06 


3/12 


8/12 


0/12 


7.4 


7.1 


NA 


26.120 


1.39 


1.06 


1/7 


5/7 


0/7 


6.8 


6.0 


NA 


27.120 


1.34 


1.05 


3/13 


0/13 


13/13 


5.5 


5.8 


NA 


28.120 


1.28 


1.03 


1/5 


5/5 


0/5 


6.9 


3.9 


NA 


29.120 


1.27 


1.03 


0/12 


10/12 


0/12 


6.0 


5.0 


NA 


30.120 


1.23 


1.02 


0/9 


9/9 


0/9 


7.1 


6.8 


NA 


31.120 


1.23 


1.02 


3/15 


8/15 


0/15 


6.6 


2.4 


NA 


32.120 


1.22 


1.02 


0/10 


10/10 


0/10 


7.1 


4.1 


NA 


33.120 


1.16 


1.00 


0/8 


6/8 


0/8 


7.4 


5.3 


NA 


10.80 


1.14 


0.99 


2/8 


6/8 


0/8 


7.8 


6.5 


3/30 


11.80 


1.12 


0.99 


0/12 


12/12 


0/12 


7.4 


6.2 


0/30 


34.120 


1.12 


0.99 


1/5 


4/5 


0/5 


6.1 


5.8 


NA 


35.120 


1.12 


0.99 


1/7 


6/7 


0/7 


5.8 


1.8 


NA 


36.120 


1.12 


0.99 


0/8 


6/8 


0/8 


5.8 


2.2 


NA 


37.120 


1.10 


0.98 


0/12 


5/12 


7/12 


5.2 


4.4 


NA 


38.120 


1.09 


0.98 


3/9 


4/9 


0/9 


7.3 


3.6 


NA 


39.120 


1.08 


0.98 


4/9 


2/9 


7/9 


8.4 


0.7 


NA 


12.80 


1.01 


0.95 


0/9 


9/9 


0/9 


8.1 


3.9 


0/34 


40.120 


1.01 


0.95 


0/6 


2/6 


0/6 


7.1 


6.1 


NA 


41.120 


1.01 


0.95 


0/9 


9/9 


0/9 


8.5 


4.7 


NA 



Note - The first column denotes the our id numbers, ordered by mass within their respective simulation box. The number following the decimal 
point corresponds to the box size in comoving fe^'Mpc. The second column lists the cluster virial mass in units of IO'^'/i^'Mq. The third column 
lists the cluster virial radius in h^'Mpc. The fourth column, A/grp/Motai, lists the fraction of {Mi„ > 1O"^/i"'M0) galaxies that were accreted 
in a group-size halos with M > IO'^/z^'Mq. The fifth and sixth columns list the fraction of galaxies that were accreted as the only object in 
their halo (No/Nmai) and with five or more companions in their halo (A'>5/A'iotai), respectively. Columns seven and eight show the approximate 
lookback times to the host halos accretion of half of its mass and half of its surviving substructure respectively. Finally, the last column is 
the same as the fourth column, except now we track the fraction of galaxies accreted in group-mass halos using our high-resolution sample of 
M,„ > 1O"/!"'M0 galaxies. 



